Shock Structures and Velocity Fluctuations in the 
Noisy Burgers and KdV-Burgers Equations 

Hidetsugu Sakaguchi 
Department of Applied Science for Electronics and Materials, 
Interdisciplinary Graduate School of Engineering Sciences, 
Kyushu University, Kasuga, 816-8580, Japan 

February 2, 2008 



abstract 

Statistical properties of the noisy Burgers and KdV-Burgers equations are nu- 
merically studied. It is found that shock-like structures appear in the time- 
averaged patterns for the case of stepwise fixed boundary conditions. Our re- 
sults show that the shock structure for the noisy KdV-Burgers equation has an 
oscillating tail, even for the time averaged pattern. Also, we find that the width 
of the shock and the intensity of the velocity fiuctuations in the shock region 
increase with system size. 



1 Introduction 

The Burgers equation and the KdV-Burgers equation are well-known model 
equations used in the study of shock waves in fluids and plasmas. 1 , 2 Spatial 
integration of the noisy Burgers equation yieldsthe Kardar-Parisi-Zhang equa- 
tion. The KPZ equation has been intensively studied in the context of growing 
rough interfaces. !3] The Kuramoto-Sivashinsky equation is one of the simplest 
partial differential equations that exhibit spatiotemporal chaos. There is a con- 
jecture that the statistical properties of the Kuramoto-Sivashinsky equation on 
large scales are closely related to those of the noisy Burgers equation. |2 El In 
previous studies, we have found shock structures in the time averaged patterns 
of the Kuramoto-Sivashinsky equation with fixed boundary conditions, where 
the boundary values are set as u(0) = —Uq, u{L) = C/q-EICI In this paper, we 
study some statistical properties of the noisy KdV-Burgers equation and the 
noisy Burgers equation with stepwise fixed boundary conditions. 
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2 Time-averaged shock structure and velocity 
fluctuations 



The noisy KdV-Burgers equation in one- dimension has the form 

Ut = Uxx + dUxxx + UUx+^xix,t), (1) 

where u(x,t) is interpreted as a velocity variable on the interval x G [0, L], 
and the subscripts denote differentiation. The noise ^{x,t) is assumed to be 
Gaussian white noise satisfying 

{^{x, t)) = 0, {^{x, mx,' , t')) = 2D6{x - x')6{t - t'). 

If the parameter d for the dispersion term in Eq. (1) is equal to zero, this 
equation reduces to the noisy Burgers equation. The quantity u{x)dx is 
conserved during the time evolution. 

We consider firstly the case of periodic boundary conditions. The Fourier 
amplitude Uk„{t) = {1/\/L) u{x) exp{—iknx)dx = 2im/L) satisfies the 
Langevin equation 

= -kluk„ + iknCknit) + {-idkluk„ + -h= X^(*fcm)wfe„-fe„.Wfe^}, (2) 

^ m 

where Cfe„(^) = (1/v^) J^^ ^{x) exp{—iknx)dx is the Fourier amplitude for the 
noise £,{x,t), which satisfies {^kj = 0, (C/c„ (i)C/c„. (*')) = '2DSn+m,oS{t - t'). 
Keeping only the first and second terms on the right-hand side of Eq. (2), we 
obtain the simplest linear Langevin equation. Both the dispersion and the 
nonlinear terms in the brackets on the right-hand side in Eq. (2) conserve 
the energy integral \u{x)\'^dx = The probability distribution 

P({?ife„}) cx cxp(— |Mfe„p/2D) is therefore an equilibrium distribution for 
the corresponding Fokker-Planck equation, 



dP ^ d 



dt ^ duk„ 



(^-f^luk^ - idkluk^ + -^^ikmUkn-k,„Uk^ P + 2Dkl 



dP 



(3) 

In real space, the equilibrium distribution can be rewritten P{{u{x)}) oc 

exp(— u{x)'^dx/2D). We have confirmed this thermal equilibrium distribu- 
tion in numerical simulation, using the Heun method (the second order Runge- 
Kutta method) with temporal and spatial step sizes At = 0.00025 and Ax = 
1/4. The parameter values are chosen to be L = 400, d = 2 and _D = 1/4. The 
numerical simulation was performed up to time tf = 101250, and the statisti- 
cal averages were calculated as the long time averages between t = 1250 and 
t = 101250. The equilibrium distribution for the spatially discretized system is 
written P{{ui}) oc exp(— u^Ax/2D), where Ui represents u{x) at a; = iAx. 
The average value of u{x) is zero, and therefore the variance of the velocity is 
{u{xy) = Dj Ax for any x. Figure 1(a) displays the profile of (^^(a;))^/^. The 
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Figure 1: (a) Time-averaged profile of {u^)^/"^ for the noisy KdV-Burgers equa- 
tion with d = 2 under periodic boundary conditions, (b) The velocity distribu- 
tion P{u) at a; = L/2. 



average value of v? is almost 1 for any x. Figure 1(b) displays the velocity 
distribution P{u{L/2)) dX x = L/2 and it is compared with the Gaussian distri- 
bution with variance 1. No difference can be discerned in this plot. The thermal 
equilibrium distribution was also obtained for the noisy Burgers equation with 
d = 0. 

Next, we performed numerical simulations with stepwise fixed boundary con- 
ditions. The boundary values we used arc u(0) = —Uq and u{L) = Uq for 
L = 50. The initial condition is u{x) = Uq ta.iah{Uo{x — L/2)/2). During the 
simulation, at regular intervals, we confirmed that the condition udx = was 
satisfied. Figure 2 displays the numerical results for the noisy Burgers equa- 
tion with Uo = 1.5. The solid curve in Fig. 2(a) represents the time-averaged 
profile {u{x)), and it is compared with the curve of 1.5tanh0.71(x — L/2). A 
shock structure appears in the time-averaged pattern. We have checked that 
the time-averaged pattern depends only weakly on the total time tf. If the 
noise term is absent, i.e., f = 0, the shock solution to the Burgers equation is 
u{x) = C/q tanh{C/o/2(.T — L/2)}. The width of the time-averaged shock struc- 
ture in Fig. 2(a) is almost the same as that of the deterministic equation. Figure 
2(b) displays a profile of the standard deviation (i5w^)^/^ = {{u{x)— {u{x)))'^)^^'^ . 
In the flat region, distant from the shock region, the standard deviation of the 
velocity distribution is almost 1, but, it increases in the shock region. The 
variance of the velocity distribution is interpreted as an effective temperature. 
If the position of the shock is assumed to fluctuate due to the noise term, 
the profile of u would fluctuate as u{x) ~ Uotanh{Uo/2{x — L/2 — Sx)} ~ 
?7otanh{[/o/2(a;-L/2)}-C/^fe/[2cosh^{C/o/2(a;-L/2)}]. Then, the deviation of 
((5w2)i/2 from 1 can be approximated as {U§/{2cosh^{Ua/2{x - L/2))}{Sx'^y/'^. 
The numerical result for the standard deviation of the velocity distribution is 
fit well by the form 1 + 0.14/ cosh2{0.75(a; - L/2)}, as shown in Fig. 2(b). The 
slight deviation of the time averaged profile from Uo tanh{t/o/2(.x — L/2)} may 
also be interpreted as a blurring effect resulting from the fluctuations of the 
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Figure 2: (a) Time-averaged pattern of u{x) (solid curve) for the noisy Burgers 
equation under the stepwise fixed boundary conditions m(0) — —1.5 and u{L) = 
1.5 with L = 50. The dashed curve represents 1.5 tanh(0.71(a; — L/2)). The 
two curves are almost indiscernable. (b) Time-averaged profile of (Su^)^^'^ = 
{{u^ix)) - (w(a;))2)i/2 for L = 50. The dashed curve is 1 + 0.14/ cosh2{0.75(a; - 
L/2)} (c) Time averaged pattern of u{x) (solid curve) and 1.5 tanh(0.475(a; — 
L/2)) (dashed curve) for L = 400. (d) Time averaged profile of (Su^)^^^ = 
iiu^x)) - (m(x))2)1/2 for L = 400. 



shock position. (This type of deviation of the time-averaged profile from the 
deterministic shock is observed for sharp shocks, for which Uq is sufficiently 
large. Such a deviation was not clearly observed in a previous investigation in 
the case of smooth shocks. 0) 

The profile of the time-averaged pattern and the standard deviation does 
not depend on the the total simulation time tf. However, it depends strongly 
on the system size, L. Figure 2(c) displays the time-averaged profile (u(x)) 
for L = 400, and it is approximated as 1.5tanh0.475(x — L/2). The width of 
the shock structure is larger for this larger size system. Figure 2(d) displays 
the standard deviation of the velocity distribution for L = 400. We see that 
the increase of the effective temperature in the shock region is larger than in 
the case of i = 50. The time-averaged shock structure can be approximated as 
1.5 ta,nh.{K{x — L/2)}, where k is the inverse of the shock width. The system size 
dependence of k and the excess part of the standard deviation of the velocity 
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Figure 3: (a) System size dependence of k, where the time averaged shock 
structm-e is approximated as 1.5 tanh{K(a; — L/2)}. (b) System size dependence 

1 /2 

of the excess part 6Te of the standard deviation of the velocity distribution. 
The dashed hne corresponds to a hne of L^/^. 

distribution StI/'^ = (5u2)i/2 _ at a; L/2 are displayed in Fig. 3. The 
quantity seems to increase as 5Te ^ L^/^. We do not yet understand the system 
size dependence well for this noisy Burgers equation. A simple interpretation 
is as follows: the value of u fluctuates randomly around ±[/o in the region 
distant from the shock position, therefore, the integration / u{x)dx between 
and L except for the shock region takes a random value of 0(L^/^). The 
total value of jj^ u{x)dx is conserved to be in our model. To compensate the 
total fluctuations of u{x) in the surrounding region, the position of the shock 
moves necessarily by — 0(L^/^). Shock fluctuations have been studied in an 
asymmetric simple-exclusion process with a blockage site, which is one of the 
simplest of the driven diffusive lattice gas madels. In that model, the amplitude 
of the fluctuations of the shock position increases as L^/^ or L^/^, depending 
on the system parameters. ;8 This behavior may be related to our results. (In 
contrast, for an asymmetric simple-exclusion process with open boundaries, the 
shock position exhibits a random walk.*^) 

Figures 4(a) and (b) display the time-averaged profiles in the cases L — 
and 400 for the noisy KdV-Burgers equation with d = 2. We see that the time- 
averaged profile exhibits damping oscillations for x > L/2, and it is asymmetric 
about the center of the shock structure. This is characteristic of the KdV- 
Burgers equation with a dispersion term. The dashed curve is the numerically 
obtained shock solution for the KdV-Burgers equation (1) without the noise 
term for d = 2. It is seen that this curve nearly coincides with the time- averaged 
structure for L = 50. However, the width of the time-averaged shock structure 
is longer for L = 400, and the fit to the deterministic shock is worse for the 
larger system. This may also be due to the fluctuations of the shock positions. 
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Figure 4: (a) Time averaged-pattern of u{x) for the noisy KdV-Burgers equation 
with d = 2 under the stepwise fixed boundary conditions u{Q) = —1.5 and 
u{L) = 1.5 with L = 50. The dashed curve is the stationary solution for the 
KdV-Burgers equation with d = 2 without the noise term, (b) Time-averaged 
pattern of u{x) for the noisy KdV-Burgers equation with L = 400. The dashed 
curve is the stationary solution for the deterministic KdV-Burgers equation, 
(c) Time-averaged profile of {du^)^/"^ = {{u'^{x)) — (u(a;))^)^/^ around the shock 
positions for L = 50 (dashed curve) and L = 400 (solid curve) . (d) The velocity 
distribution P('u) at a; = L/2-1, L/2, L/2 + 1, L/2 + 2 and L/2-F3 for L = 400. 
The numbers in the figure denote the distance, x — L/2, from the center. 

Figure 4(c) displays the standard deviation of the velocity distribution {Su"^)^^"^ 
for L = 50 (dashed curve) and 400 (solid curve). It is seen that here too, the 
effective temperature increases in the shock region. The profile of the effective 
temperature is also asymmetric around the shock center. The increase of the 
effective temperature for L = 400 is clearly larger than for L = 50. Figure 4(d) 
displays the velocity distribution P{u{x)) at a; = L/2 — 1, L/2, L/2 + 1, L/2 + 
2 and L/2 + 3 for i = 400. The velocity distribution approaches Gaussian 
distribution with variance 1 as \x — L/2\ is increased. The variance of the 
velocity distribution is larger than 1 and the distribution becomes asymmetric 
near the shock region. 
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Figure 5: (a) Energy dissipation rate (w^) for the KdV-Burgers equation with 
d = 2 and system size L = 400 under periodic boundary conditions, (b) En- 
ergy dissipation rate (sohd curve) and the contribution from the velocity 
fluctuations {{Su)-^) (dashed curve) for the KdV-Burgers equation with d = 2 
under the stepwise fixed boundary conditions u{0) = —1.5 and u{L) = 1.5 for 
L = 400. 

3 Energy dissipation rate 

The nonequilibrium states discussed above appear as a result of the stepwise 
fixed boundary conditions. Energy is injected into the system through the 
boundary conditions, and energy dissipation occurs in the region of the shock. 
From Eq. (1), the time evolution of the energy density can be written as 



1 du d{uux) 
2~dr ^ dx 



dx 2 dx 



-^-ul + u^.. (4) 



If D = and Uq 0, the third term on the right-hand side is related to the 
energy injection at the boundaries, and the fourth term represents the energy 
dissipation. The energy injection from the boimdary conditions is 2/3J7g. The 
total energy dissipation is expressed as Q = Jq u'^dx, and it is equal to 2/3[/q. 
(This can be explicitly shown by the integration of = {dx{Uota.nhUo{x — 
L/2)/2})^ for the case d ~ 0.) If D 7^ and periodic boundary conditions 
are employed, the total energy dissipation rate, Q = u^dx, is equal to the 

quantity Q' — u^^dx, which is interpreted as the heat flow from the heat 
reservoir expressed by the noise term. The total entropy production is therefore 
zero in the thermal equilibrium state. We have calculated the time average 
{u^(x)) in our simulation model with Ax =1/4 and d — 2. Figure 5(a) displays 
the profile of (u^(a;)) for the case of periodic boundary conditions. In this case, 
the energy dissipation rate is nearly 32. Also, it can be shown that the average 
{u{x,t)£,x) = 2D/{Ax)^ ~ 32 in the spatially discretizcd system with D = 1/4 
and Aa; = 1/4. That is, we have found numerically that the average of the 
entropy production rate is zero for any x. The probability distribution of 
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is well approximated as a Gaussian distribution. The solid curve in Fig 5(b) 
displays the profile of {u'^{x)) for the noisy KdV-Burgers equation in the case 
of c? = 2 with the fixed boundary conditions u(0) = —1.5,u{L) = 1.5 for L = 
400. Note that the energy dissipation rate is larger than the equilibrium value 
32 in the shock region. The energy dissipation term is written as a 

sum of two components: u'^{x) = + {Su{x))'^. The dashed curve in 

Fig. 5(b) represents the contribution by {(Su)"^) from the velocity fluctuations. 
The rate of energy dissipation resulting from the velocity fluctuations seems to 
be larger than that resulting from the time-averaged pattern in this simulation. 
We calculated the energy dissipation rate resulting from the velocity fluctuations 
also for L = 50. Clear difference by the system size was not found in this 
energy dissipation rate resulting from the velocity fluctuations in contrast to 
the velocity fluctuations shown in Fig. 4(c). 

4 Summary and discussion 

We have performed numerical simulation of the noisy Burgers and KdV-Burgers 
equations. We have found that a thermal equilibrium state is realized for 
periodic boundary conditions and a shock structure appears as a nonequilib- 
rium state for stepwise fixed boundary conditions. The statistical properties 
of nonequilibrium states have been studied by many researchers. In particu- 
lar, the statistical properties of nonequilibrium liquids have been intensively 
studied using molecular dynamics simulations. |10| Nonequilibrium states, in- 
cluding shocks, have been studied using driven lattice gas models, most notably 
those with asymmetric exclusion processes. [TTl 15] We have studied the statisti- 
cal properties of the nonequilibrium states created by the boundary conditions 
with our Langevin-type models. Increases in both the effective temperature and 
the energy dissipation rate caused by the velocity fluctuations are observed in 
the shock region. This increase of the effective temperature may be due to the 
fluctuations of the shock position. The intensity of the fluctuations depends 
strongly on the system size. The increases of the energy dissipation rate re- 
sulting from the velocity fluctuation does not so depend on the system size. 
The mechanisms of the increases of the two statistical quantities in the shock 
region may be different. The time- averaged shock structure in our model may 
be closely related to the macroscopic shock structure found in the asymmetric 
simple exclusion model. However, a detailed investigation of their relation is 
left to a future study. 
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